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ABSTRACT 

A summary of QCD results obtained at e"^e~ experiments in recent years is given. 
Emphasis is put on basic QCD tests and measurements with event shapes. 
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1 Introduction 

Experimental studies of strong interaction physics at e^e^ machines aiming for basic 
tests of Quantum Chromodynamics, QCD, enfold the determination of the coupling 
and the quark masses, the parameters of the Lagrangian , the verification of the 
basic vertices of QCD, and tests of the loop corrections interrelated with asymptotic 
freedom and confinement. Moreover the experimentalists are faced with a wealth of 
strong interaction physics phenomena. 

e^e~ annihilation to hadrons provides the simplest strongly interacting ini- 
tial state, a quark anti-quark pair. Actual measurements span from the $ resonance 
at Da$ne to ^/s ~ 209 GeV at Lep II. This talk reviews some basic measurements 
performed by e^e~ experiments in recent years. 

Results still come up from the measurements at the Z where the cross- 
section is large and background free precision measurements can be performed. At 
higher energy the measurements are complicated by frequent initial state photon ra- 
diation, leading to "Z return" events, and boson {W^W~ and ZZ) pair production. 
Still the LEP data in combination with previous low energy results allow detailed 
tests of the energy evolution of hadronic properties. 

Off resonance data from the B factories are important for studies of charm 
final states as, in this case, the charmed hadrons are not influenced by B hadron 
decay products. Moreover radiative Onia decays so far provide the only lucid source 
of purely gluonic final states. 

2 Energy Dependence of the Hadronic Multiplicity 

The dependence of the charged hadron multiplicity on the CMS energy up to the 
highest Lep II energies is shown in Fig. |3 The basic process governing the multi- 
plicity increase is gluon Bremsstrahlung off the initial gg pair (oc Cp-ois] Cp = 4/3). 
The measurements are well described by parton shower (PS) fragmentation models, 
as well as by Modified Leading Log Approximation (MLLA) and 3NL0 calcula- 
tions pp. The calculations intrinsically predict the number of final state gluons and 
quarks, the connection to the hadronic multiplicity being made by Local Parton 
Hadron Duality (LPHD) assuming the number of hadrons to be proportional to the 
number of partons. The proportionality factor is fitted to data, leaving only the 
energy increase predictable. 

An important common aspect of all descriptions is the suppression of soft 
gluon emission due to gluon coherence. Impressive evidence for this phenomenon 
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Figure 2: CMS energy dependence of 
the maximum of the ^ distribution, ^* . 



comes from the "hump-backed" plateau shown in Fig. [T] depicting the suppression of 
low energy hadrons. The data over a wide energy range are compared to the MLLA 
"limiting spectrum" prediction. The increase of the maximum value, of these 
spectra is strongly reduced compared to phase space expectation (Fig. 

Independent insight on gluon coherence stems from the multiplicity dif- 
ference, 6bi, between light / = uds and b quark events. In b events the collinear 
divergence of gluon radiation is regulated by the b mass. Radiation is suppressed 
within a "dead cone" of opening angle G oc mb/E, i.e. decreases towards high 
b energies. The corresponding increase in phase space for radiation and gluon co- 
herence lead to the expectation of constant 6bi- This expectation has been recently 
verified by Lep II measurements j2]. By contrast, incoherent radiation would lead 
to a strong decrease with energy: 6bi ~ 2{Nb) +Nii{{l— {xB))'^s) — Nqq{s), with {xb) 
the average energy fraction and (Nb) the average decay multiplicity of B hadrons. 

Analogously to the quark case the multiplicity in gluon jets is governed 
by gluon Bremsstrahlung off gluons (oc Cac^s)- Due to the higher colour factor 
Ca = 3, a multiplicity increase by the factor Ca/Cf = 2.25 is expected at high 
energy. Experimental results obtained mainly in e~^e~ 3 jet events were found in 
the range 1 ... 1.5. The small ratio is now understood considering gluon coherence 
(leading to p± like scales), finite-energy/non-perturbative (n.p.) effects and biases 
due to jet selection. 

Evidence for n.p. effects comes from the rapidity spectrum of hadrons 
in "leading" gluon jets recoiling with respect to a bb system of comparably small 
inter quark angle ~ 90° j^. For small rapidity, i.e. for particles produced first 
in the hadronisation, the production rate is about twice as high as for quark jets 
of corresponding energy. At high rapidity however, more hadrons are produced in 
quark jets. This is partly due to leading particle effects - quarks are "valence" 
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Figure 3: Charged hadron multiplicity of 
three jet events as function of the inter 
jet angles. 
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Figure 4: Charged multiplicity for qq and 
gg colour singlet states as function of the 
CMS energy. 



particles of hadrons. Moreover, as Bremsstrahlung is less effective in quark jets, 
more energy remains available and allows high energy hadrons to be formed. The 
multiplicity of 3 jet events is predicted as: 

Nggg = Ngg{Sqg, Vent) + ^Ngg{plg) , (l) 

where Ngg is the multiplicity observed in e~^e~ annihilation, Ngg is the multiplicity 
for a colour singlet gg pair. The choice of scales Sgg and p^^ is ambiguous. It 
depends on the part of phase space assigned to the quark and gluon jets in a 3 
jet event, respectively. The ycut dependence enters solely due to a bias induced by 
the jet assignment. For a 3 jet event resolved at a selected it is prohibited to 
have further gluons emitted at higher y^^^ > y^^^ as these would lead to y^^^ = y^^^. 
Consequently a negative bias on the multiplicity is implied. 

The energy slopes of the multiplicity of gluons and quarks are related by: 



dNgg 



dU 



L'=L+c„ 



I C - ^) i'^"-*^) — 
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The solution of this differential equation implies a const, of integration suited to 
absorb n.p. qg differences. The complete ansatz describes the observed multiplicity 
for 3 jet events of general (Fig. IHl) and symmetric topology [0]. The energy scales are 
determined from the 3 jet topology and the three possible assignments of the gluon 
jet are considered by appropriate weighting with the 3 jet matrix element. Eq. is 
easily solved for Ngg. The observed increase with energy is about twice that of the 
quark case (Fig. and can be immediately recognized as due to the higher colour 
charge of the gluon. The quark curve is a MLLA fit to the data, the gluon curve is a 
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prediction deduced from the qq data using Eq. (j21). The n.p. offset was fixed by the 
Cleo point at ^/s ~ 10 GeV measured from xb 99 decays [7j. The prediction 
agrees with the 3 jet data from Tristan and Lep jSj. From symmetric 3 jet events 
Delphi measured the colour factor ratio: 

^ = 2.221 ± 0.032,tat. ± 0.047exp. ± 0.058hadr. ± 0.075theo. (3) 
Up 

3 Heavy Quark Fragmentation 

Despite their complex decay patterns heavy hadrons can be reconstructed due to 
their long lifetime. Since heavy quarks are rarely produced in fragmentation it is 
possible to determine even the fragmentation function of primary B and D hadrons. 

The expendable micro-vertex detectors of the Lep/Slc experiments allow 
inclusive B reconstruction with high purity and efficiency at a (9(10%) energy res- 
olution. The new measurements of the b fragmentation function ^ (Fig. agree 
with each other. Note however, that the data points are correlated due to the limited 
energy resolution and the unfolding applied in the analysis. The precise average of 
the scaled energy fraction {xb) = 0.715 ± 0.003 of weakly decaying 5's is slightly 
increased compared to previous e+e^ results. This may be due to the inclusion 
of B baryons in the sample. The measurements now allow to distinguish between 
different fragmentation models. Consistency is only obtained with the Lund and 
the similar Bowler model while the data disagree with the widely used Peterson 
model. Very precise data on c fragmentation also became recently available 0. The 
Belle measurement is based on exclusively reconstructed D*^ mesons while Cleo 
also uses and decays. To avoid contamination from B decays at low xd 
off resonance data {C ~ 3fb~^) are used. For xd > 0.5, where B feed through is 
kinematically excluded, also the huge T(4S') data sets (£ ~ 6 — 25/6^^) are used. 




Figure 5: B hadron fragmentation func- 
tion as measured by Sld and Lep. 
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Figure 6: Fragmentation function of D* 
mesons as measured by Cleo. 
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The momentum resolution is excellent due to the low CMS energy. The Belle 
measurement gives {x^*) = 0.612 ± O.OOSstat ± O.OO^syst and strongly favors the 
Bowler model. Cleo observes differences in the spectra of D*^, and D'^, which 
are likely due to differing influence of D resonance decays. Such effects have to be 
considered when interpreting the precise results [Hj. 

The application of the e^e~ heavy quark results in the description of other 
processes requires a proper extraction of the n.p. fragmentation function. The e+e~ 
analyzers now perform this extraction fLU\. According to the QCD factorisation 
theorem the measured function is a convolution of a pert, piece, describing hard 
and soft gluon radiation, and the genuine n.p. fragmentation function. After Mellin 
transformation the convolution turns into a simple product of transformed contri- 
butions. This equation can be easily solved for the n.p. part and transformed back. 
The resulting n.p. B fragmentation function as evaluated from two Lep measure- 
ments agrees well ^Uj • It is strongly peaked around ~ 0.9 showing that most 
of the b quark energy is conserved in the B hadron. This function can be directly 
applied in theoretical predictions provided similar assumptions (like renormalisa- 
tion scheme and factorisation scale) are used in the application as in the e~^e~ pert, 
calculation. 

The Sld experiment measured the double inclusive B fragmentation func- 
tion JT]. This function depends on the energies of and the angle, between the B 
particles. The ratios Gmn of the double to the two corresponding single Mellin mo- 
ments were determined. These ratios are advantageous for two reasons: The collinear 
singularity is canceled by the finite b mass and large log. terms oc In Ecm/ cancel to 
all orders [T^. The prediction for Gmn{^) agrees well with the measurement |1H I12j . 
As the n.p. fragmentation function enters twice in Gmn Sld interprets their mea- 
surement as a basic test of the QCD factorisation theorem. The theoretical work 
[12] suggests Gmn for a precise extraction of as- 

4 Measuring mb(Mz) 

Besides the coupling the quark masses are the only free parameters of the QCD 
Lagrangian. The role of the masses is similar to that of the coupling. Besides the 
pole mass, Mg, therefore the renormalised mass, mq{Q), is defined, customarily in 
the MS scheme. Quarks are bound inside hadrons, therefore their masses can be 
assessed only via dynamical relations. At high energy it is possible to determine the 
mass of the b quark from the reduction of gluon Bremsstrahlung in b compared to 
light quark events. This reduction is oc ml/p^ ^ ml/{s ■ Ucut)- The llycut factor 
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increases the effect by one order of magnitude compared to the directly expected 
m1/s dependence. Several coinciding predictions exist (see [T^). 

The mb{Mz) measurements performed by most of the Lep/Slc experi- 
ments jlSI proceed as follows: Jet rates or event shapes are measured for h and light 
quark (or all) events. Heavy quark tagging is performed using impact parameter 
based techniques. Detector and hadronisation corrections partly cancel in the ratio 
i?3 of h and light distributions leading to controllable uncertainties despite the weak 
effect ~ 3% at = Mz- 

Improved understanding of the hadronisation cor- 
rection has been gained lately jEj. Apparent 
differences between models were traced to an im- 
perfect prescription of the xb distribution. The 
new Delphi analysis avoids badly described re- 
gions henceforth leaving the associated uncer- 
tainty unimportant. A toy simulation connect- 
ing the different h quark masses and effects in 
Pythia with each other showes for the first time 
the expected stability of the hadronisation cor- 
rection against parameter changes. However, 
the uncertainty of the h mass in the model, which 
could be identified with the pole mass, now ap- 
peared as important additional uncertainty. 
The expected oc y^^^ dependence of -R3 is seen [12] and the data are con- 
sistent with the NLO expectation for mb{Mz) ~ 3 GeV. This small value already 
indicates the expected "running" of the mass seen in Fig.[7| Here the average mfc(mb) 
determined from bound states is compared to the mh{Mz) results and their average: 
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Figure 7: E dependence of m^^^ 
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For the averaging, theoretical and hadronisation uncertainties were taken fully cor- 
related. The stat. error of Opal determined from a correlated average has been 
scaled according to PDG rules with \Jx^ l^df of this average. Improved measure- 
ments of mi){mh) from moments of the lepton spectra of B decays or using lattice 
calculations are expected fH]- The fall off of the mass with energy is given by the 
Renormalisation Group Equation: 



d\nm{Q^ 
ainQ2 



-7m(as) 



7n 



(5) 
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where 7^ is the mass anomalous dimension, sensitive to loop corrections to the 
quark propagator. The succesful description of the b quark mass running implies a 
basic test of the quantum loop structure of QCD. 

5 Four Jet Events and Colour Factors 

In 4 jet events all basic QCD vertices except the suppressed quartic gluon vertex 
occur as building blocks already at tree level. At tree level the 4 jet cross-section 
can be decomposed into terms proportional to the colour factor products Cp, CfCa 
and CpUfTp. This possibility is due to the general helicity structure of the process 
and allows measurements of the colour factor ratios Ca/Cf and nfTp/Cp from the 
4 jet angular distributions. This is equivalent to a determination of the gauge group 
underlying the strong interaction. Recently NLO predictions for the cross-section 
became available (for references see ^Hl) improving the data description. At the 
cost of a leading order determination it is thus possible to determine directly the 
colour factors Ca and Cj?, as well as in combination with measurements of the 3 
jet rate ^Hl- The current averages of the colour factors are: 

= 1.34 ± 0.23 Ci. = 2.97 ±0.5 (6) 

The results for the colour factor ratios are included in Fig. A drawback of all 
measurements is the residual dependence on a necessary hadronisation correction. 
Here QCD is presumed uncared if MC models or QCD power corrections are used. 

6 Us from Event Shapes 

The measurement of with event shape observables probes the amount of gluon 
radiation in the hadronic final state oc Cp ■ ots- It is insensitive to the underlying 
electroweak physics. Event shape observables are required to be insensitive to soft 
and coUinear gluon emission as well as to hadronisation effects. A typical observable 
is t = 1 — Thrust. In the following an arbitrary observable is denoted by where 
low values indicate small and large values {y ~ 0.3) imply hard gluon radiation. 

The Us measurement proceed as follows: The experimental distributions 
are determined and corrected for cuts, limited detector acceptance and resolution. 
At high energy Z return events and two boson production need special regard. 

Two different types of predictions are available based on NLO matrix 
element (ME) and Next to LLA (NLLA) calculations. The former applies to hard 
radiation, the latter to the two jet regime (small y) where multiple gluon emission 
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needs to be resummed. Both calculations can be matched in order to obtain a 
prediction valid for a wide region of phase space. 

When comparing theory to data it is important to correct the influence 
of the n.p. hadronisation. The classical correction is given by the ratio of the y 
distribution on parton (that is just before hadrons are formed) to hadron level in a 
Monte Carlo model. Alternatively so called Power corrections (PC) are used. 

6.1 a, from NLO Matrix Elements 
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Figure 9: as{Mz) as determined from 
18 observables from Delphi data using 
C(a^) theory. 



The NLO prediction for the distribution of a 3 jet like observable is: 

D{y) = -^=A{y)cysil^) + [Aiy)Po'^^x^ + B{y)]cy^^{lx) a = — (7) 
a ay vr 

where A{ii) and Biy) are pert, calculated coefficients. For small y this prediction 
receives important higher order corrections cx \v? y. Uncalculated higher oder 
terms lead to the dependence on the unknown (observable dependent) renormalisa- 
tion scale {x^ = fi'^/s). The term changes the slope of D{y), typically falling 
off exponentially at not too small y. Satisfactory agreement {x^/Ndj ~ 1) of the 
prediction is only obtained with optimised x^. Taking x^ = 1, that is the pure MS 
prediction, will in general lead to a dependence of ag on the chosen fit range. This 
fact was long known for small y but has also been verified at high y with the precise 
Z data In fact this behavior must persist at high y as a consequence of the 
normalisation of D{y) to the total cross-section (see Eq. 
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The optimised scales for observables with large positive (negative) sec- 
ond order prediction tend to be small (high). Currently it is little understood 
why the experimentally optimised scales correlate with the theoretically motivated 
FAC/PMS/ECH scale choices (Fig.|Hl [HI). The results obtained using optimised 
scales only show a scatter of about 2%. For the = 1 results, however, due to the 
aforementioned discrepancies the scatter is bigger ~ 10% for the same y intervals 
in the fits. In this case an additional theory uncertainty is required in order to 
assure consistency of the results obtained from different observables. 

6.2 Us from Matched 0{a1) ® NLLA Predictions 

The NLLA prediction has the general form: 
ry 

R{y)= j D{y)dy = exp{Lgi{asL) + giiasL)} L = -\ny . (8) 

The function R{y), however, does not obey the phase space boundary conditions 
RiVmax) = 1, dyR{ymax) = 0. This imperfection is corrected by the transformation: 



L' = - In 

P 



' +1 



max J 

where Xl and p are arbitrary parameters [TH] . This transformation, though primarily 
modifying the prediction at high y also implies a % level shift in the as fit range. 

In order to obtain a prediction valid over a large range of phase space the 
NLLA and NLO ME predictions for a;^ = 1 are merged. Duplicated terms are sub- 
tracted by expanding the NLLA prediction in a power series in a^. This subtraction 
is ambiguous as it can be performed for i? or Ini? (see Eq. ((HI)), respectively. 

In Fig. ^1 a comparison of the T distributions measured between the Z 
and 206 GeV and the matched prediction is shown. The data are well described over 
the full range of energy and T. Only for the precise Z data a slight difference in the 
slope is indicated. This difference is weaker for most other observables. 

The Lep QCD Working Group attempts to combine the as measurements 
from the high energy ( low statistics ) Lep II measurements. The 0{a1) © NLLA 
prediction is use for the observables 1 — T, M/^, C, Bt, and ?/3. The theory 
implementations and hadronisation corrections have been compared, the treatment 
of correlations and the theory uncertainty was discussed. It was seen that often badly 
described observables dominate in the combination and a downward bias in as was 
observed. This bias is due to the proportionality of the theory uncertainty to a^. 
To avoid these faults a novel definition of the theory uncertainty, the uncertainty 
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band method, was introduced: the matched reference prediction with a fixed initial 
as{Mz) was varied by changing the renormahsation scale (1/2 < < 2), not 
applying or applying the NLLA phase space condition with 2/3 < xl < ?>l2,p = 1;2 
and changing from Ini? to i? matching. The envelope of the variations defines 
the uncertainty band which in turn determines the uncertainty of the reference a^. 
Using this fully correlated uncertainty as well as the experimental and hadronisation 
uncertainties the Lep average is then calculated. The procedure is iterated until 
the result is stable. The results for as{Mz) are P^: 

LEP I : 0.1200 ± 0.0002,^* ± O.OOOSe^p ± 0.0010/,ad ± 0.0048t,, (10) 
LEP II: 0.1201 ± 0.0005,tat ± O.OOlOe^p ± 0.0007^ ± 0.0045*/, (11) 

The observed energy dependence of shown in Fig. ^2 agrees well with the ex- 
pectation. Despite the strongly different statistics the overall uncertainties of the 
Lep I and II results are comparable. This is due to the dominance of the theory 
uncertainty and a decrease in the hadronisation (oc V-E) and theory uncertainties 
(oc al{E)) with energy. The total uncertainty of the order 4% is confirmed by the 
scatter of the results from the six individual observables studied. 

6.3 Power Corrections (PC) 

Assuming for the hadronisation process a longitudinal, "tube-like" phase space with 
fixed average transverse momentum, where only the rapidity range expands with 
increasing energy, i?, it is easily shown that hadronisation corrections vanish (x 1/E. 



11 



Also the replacement of the unknown analytical behavior of the strong coupling 
below an IR matching scale /i/ ~ 1 — 2 GeV by an effective average value Q;o(yU/) 
leads to a inverse power law behavior pUj : 



ao(/i/) - asifi) - ( 60 ■ log ^ + — + 26o ) a^(yu) 



For most observables the value y of the pert, calculation is just shifted by +CyV. 
Correspondingly its mean value {y) gains an additive correction. The constants Cy 
have been evaluated for NLLA resummable observables. The overall size of the PC, 
i.e. ao, is determined from a fit to data. Within the uncertainty of the calculation 
(~ 20%) universal cto (and ag) values are expected from the different observables. 

The model scetched above is able to succesfuUy describe the data on dis- 
tributions and means over a wide energy range |221- However, the obtained 
values of ao ~ 0.5 are only consistent within the quoted uncertainty of the model. 
Specifically ao is smaller for mean values than for distributions. Due to the obvious 
correlation of and (see Eq. f|T^ ) this leads to a corresponding difference in 
as{Mz)- The Jade experiment for instance reports 



means: a^iMz) = 0.1187 ± O.OOMfit ± O.OOOlsystaoo^Lh (13) 



sys-0. 0015th 

distributions: ^.(M^) = 0.1126 ± O.OOOSfit ± O.OOSTsystaoeoth (1^) 



This discrepancy persists even if hadron mass induced PC's are included pT|. 
6.4 Overview of Results from LEP 

Fig. ^1 compares some important results obtained from Lep data. All results 
are compatible with the current PDG average. The indirect measurements from the 
Z shape and Rr are slightly above the average. The scatter of these partly highly 
correlated results obtained using NNLO theory confirms their error estimates. 

The classical ( 3 & 4 jet ) NLO direct measurements (with and Monte 
Carlo hadronisation) agree well with the PDG average. The highly correlated result 
using matched theory is slightly higher also in comparison to the pure NLLA result. 
A shift towards higher can be understood if the (9(a^)|a;^=i part of the matched 
prediction underestimates the slope of the data distribution. The results for means 
using PC's are slightly higher than the PDG average. Overall the scatter of the 
direct results confirms an error of 3 — 4% of Og typically assigned to the individual 
measurements. Improved predictions (like complete 0{al) ME's) are mandarory for 
any significant reduction of the uncertainty. 
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Figure 12: Comparison of some impor- 
tant as results from Lep. 

6.5 The Running of ag 



Figure 13: Energy evolution of (1 — T). 
The slope determines the (3 function. 



The striking strong interaction physics phenomena of asymptotic freedom and con- 
finement are intimately coupled to the energy evolution, the "running" of as- In 
QCD this is controlled by the Renormalisation Group Equation (RGE): 



das 



aing2 



The (3 function describes the loop corrections to the gluon propagator, to leading 
order as: 

It is influenced by all strongly interacting particles. A measurement of the (5 func- 
tion implies model independent limits on hypothetical strongly interacting particles. 
Such a measurement has recently been performed ^21] based on the Renormalisa- 
tion Group Invariant (or Improved) pert, theory (RGI p.t.). It applies to observables 
Roc {y)/A depending on a single energy scale which are normalised such that their 
pert, expansion starts with a^. RGI p.t. requires the observables to fulfil the RGE. 
This equation is solved as for as yielding an R dependent integration constant, A^^, 
which can be determined from a fit to the data. Ar is exactly related to ■ 
RGI p.t. is numerically equivalent to the theory in the Effective Charge (ECH) 
scheme. E.g. RGI determines A^^^ without freedom of renormalisation scheme 
choice. It is important to note that RGI p.t. implicitely resums UV divergencies. 
Power corrections can be included. 
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The RGI prediction has been fitted to the energy dependence of 7 event 
shape means |21J. The description of the data is equally good as with PC's though 
with less parameters. For RGI the power terms turned out to be compatible with 
zero (< 2% at the Z) and in contrary to the MS+PC case all observables can be 
described with a single value of as = 0.1201 ±0.0020. I.e. the energy dependence of 
mean values can be described at the 2% level without hadronisation correction and 



without renormalisation scale freedom. This implies that in the MS scheme power 
terms for mean values to a large part parameterise missing higher order terms. In 
fact it is even possible to predict the "n.p." terms using RGI [2l|. 

To leading order the RGE (Eq. (jl5|) ) specifies a straight line when plotting 
the inverse of the observable vs. the logarithm of the energy, InQ. Thus, up to a 
small correction the slope of the data directly determines /3o- In order to optimally 
determine the f3 function data on (1 — T) has been chosen spanning the energy range 
15 to 205GeV(Fig.[TnD. A linear fit yields: 

^^^^8.70 ±0.35 , (17) 

consistent with the expected QCD value 8.32. This result model independently 
excludes supersymmetric gluinos up to masses of 30-40 GeV ^Tj. Assuming QCD 
/3o or can be alternatively determined: 

/?o = 7.83 ± 0.32 n/ = 4.75 ± 0.44 . (18) 

This result is a clear measurement independent of assumptions about the /3 function 
unlike the case when using models in determinations. It is superior to results 
from Lep (/3o = 7.67 ± 1.63) HHf or world data on a, (/5o = 7.76 ± 0.44) |2S!. 

7 Summary 

Measurements performed at e^e^ machines, especially those at and above the Z 
pole over the past decade provide a vast amount of hadonic physics results. For this 
talk the selected results appeared during the last few years. 

Understanding of the dynamical behaviour of hadrons has been attained. 
Most important ingredience of the description here is the coherence of gluon radia- 
tion. Moreover it has recently been possible to verify the different coupling strengths 
of quark gluon and gluon gluon vertices directly from the hadrons created in three 
jet events. 

The instrumentation of the experiments allowed detailed measurements of 
heavy quark production. Measurements of the heavy quark fragmentation function 
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will serve as an important input to hadron colliders. The measurement of the running 
of the b quark mass tests the loop structure of the b quark propagator. 

Several types of measure- 
ments performed lead to a consistent re- 
sult for the strong coupling, though with 
a still sizable uncertainty at the few % 
level. In most cases the uncertainty is 
dominated by the theory uncertainty. For 
a sizable reduction in as error improved 
predictions are required. 

The measurements of the colour 
factors from 4 jet events, the multiplicity 
ratio obtained in gluon and quark jets 
and the quantum loops verified by the 
measurement of the /3 function strongly 

Figure 14: Caszmzr factor raUos for ^^^^^^^ §^^8^ S^^^P ^^^^^^g i^^^^" 
several groups compared with constraints actions to SU{3) / QCD (Fig. I14j). 
from various measurements. 
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